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Abstract 

The intermittency analysis of single event data (particle moments) in mul- 
tiparticle production is improved, taking into account corrections due to the 
reconstruction of history of a particle cascade. This approach is tested within 
the framework of the a-model. 



*e-inail: uf r j anikO j etta . if . uj . edu . pi, beatazSqcd . if j . edu . pi 
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1 Introduction 



The first data on possibile intermittent behaviour in multiparticle production [|l| 
came from the analysis of the single event of high multiplicity recorded by the 
JACEE collaboration |0. It was soon realized, however, that the idea may be ap- 
plied to events of any multiplicity provided that averaging of the distributions is 
performed [Q. This led to many successful experimental studies of intermittency 
[Q, and allowed to express the effect in terms of the multiparticle correlation func- 
tions 0. It should be realized, however, that the averaging procedure, apart from 
clear advantages, brings also a danger of overlooking some interesting effects if they 
are present only in a part of events produced in high-energy collisions. For example, 
the unique properties due to the presence of quark-gluon plasma in multiparticle 
production would manifest only in some events, see e.g. f^. Taking into account 
the sample of events and averaging over them destroys such an information. There- 
fore, as already discussed in §|, 0, there is a need for event-by-event analysis of 
multiparticle production processes. In this way the fluctuations of the measured 
physical quantities (e. g. factorial moments) from event to event can be observed 
and estimated, and any anomalous behaviour of them has a chance to manifest very 
clearly. Such studies should necessarily be restricted to high-multiplicity events 
because only there one may expect the statistical fluctuations to be under control. 

Such an approach to the multiparticle data analysis has been already proposed 
in [0, [|], 1^. In [H a new quantity: erraticity has been introduced to investigate the 
event-by-event fluctuations of factorial moments, and to search for their properties. 
Erraticity denotes the normalized moment of event-by-event distribution of a hori- 
zontally averaged factorial moment. It probes both types of fluctuations: horizontal 
ones connected with the spatial bin pattern and vertical ones i.e. event-by-event 
ones. 

In the event-by-event fluctations of particle moments have been investigated 
directly for the one- dimensional a— model of random cascading. Monte-Carlo simu- 
lations of the model allowed one to obtain the histograms of event-by-event distribu- 
tions of horizontally averaged particle momenta and estimate the relation between 
the intermittency parameters obtained from such a histogram, and the intermit- 
tency parameters derived after usual procedure of averaging particle moments over 
all events. The results were promising: the average value of the intermittency expo- 
nent reproduced well the value obtained by averaging particle moments over events, 
however with the tendency to underestimate the theoretical value. Furthermore, 
the dispersion of the moment distribution was inversely proportional to the length 
of a generated cascade, and even for short cascades substantially smaller than the 
average value. The latter property was of a special importance : it allowed one to 
distinguish between groups of events emerging from cascades with different charac- 
teristics. 
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In this paper we would like to improve the analysis of single event data presented 
in Taking into account corrections due to the method of recovering the history 
of the multiparticle cascade 0, we expect to reduce the discrepancy between the 
theoretical value of intermittency exponent and its value estimated from the event- 
by-event histogram |^ . Our discussion will proceed as follows. In section 2 we recall 
the definition of the intermittency exponents and the technique used to calculate 
them , 0] . In section 3 the definition of the a model will be briefly presented, and 
applied in section 4 to calculate corrections for extracting intermittency exponents 
from single event data. Section 5 is devoted to the comparison of theoretical results 
with numerical simulations. Finally in section 6 we present our conclusions. 



2 Intermittency exponents 

Consider a multiparticle production cascade distributed into M bins. At the n*'* 
stage of the cascade we measure the distribution of a particle density into M bins. 
Assume for simplicity that M = 2". We thus have 2" numbers (quantities) denoting 
the content of each bin : 

xt''> ,t = 0,l,...X-l. (1) 

To perform the event-by-event analysis one is interested in the behaviour of particle 
moments with the stage of the cascade : 

4"' = (-!"')'■ (2) 

^ i=0 

The scaling behaviour of these moments is parametrized by intermittency exponents 

_ 2"-<^^ (3) 

The task is to estimate the value of an intermittency exponent. There are two 
different ways of doing it. The first one is to calculate the average moment zjp^ for the 
whole ensemble of individual events, and from this to reconstruct the intermittency 
exponent. The second one is to calculate the exponent (pg for each event separately, 
and then to recover the average (pg. The latter approach has the advantage of being 
able to distinguish between two independent cascading processes each with different 
(j)g. This could be done by looking at the distribution of individual 0g's. In the 
former method both of these possibly independent processes would be artificially 
forced to be described by a single 'effective' (j)g. 

In the following we would like to address the question of reliably reconstructing 
the correct value of (pg from single event data. Numerical simulations in |^ showed 
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that there is an inherent discrepancy between the theoretical value and the distri- 
butions of event-by-event (pg (see Tables 1,2). The aim of this letter is to analyze 
this result and introduce a correction which improves the estimation. 

A convenient way of calculating 0^ is to make a linear fit to the points {n, log z^^^) 
( all logarithms are taken to be calculated in base 2, i.e. logx = lnx/ln2) : 

log4") =n-0, + 6. (4) 

This procedure has the advantage of cancelling out the major part of the correc- 
tion coming from the fact that we are effectively reconstructing the exponents from 
(logzjp'^'j while the true value is defined in terms of log l^z^^^^. 

However there is still one caveat to (^. Since we cannot in general separate out 
the various stages of the cascade, one reconstructs the previous stages from the last 
one by summing the x-^^'s in adjacent bins using the technique described in |I| ( 
and applied there to JACEE event ). Namely one approximates the true value 
of by : 

(n—k) (n—k) (n) 

^ j=0 

Therefore in (HI) one really uses the reconstructed moments : 



z. 



(fc) - A f.S^A 



q, reconstructed nh 

We will now use the a model of random cascading to calculate explicitly the 
difference between the true and reconstructed moments and the resulting shift of 
the (pq distribution from the theoretical value. 



3 The a model of random cascading 

In the a model of random cascading [|1| the root of the cascade — Xq*^ is taken 
to be a with probability pa and b otherwise (with probability pb = I — Pa)- One 
generates the next stages of the cascade recursively. The two bins x^^^^ and x'^^^l^ 
are obtained from x\"^ by : 

X2i~^^^ — ^ ■ a;^-"'' with probability pa, (7) 

^2i~^^^ — ^ ^ ■ ^i^^ with probability pb, (8) 

and same for x''2i^l\ The parameters a and b are taken to satisfy : 

apa + bpb = 1. (9) 
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Particle moments fulfill the relation : 

^(n) _ 2("+l)-<^^ (10) 

where intermittency exponents (f)q are equal to : 

0g = log(aV + 6%). (11) 

4 Reconstructed moments 

The reconstructed moments in the a model are related to the true ones by : 

^q,reconstructed ^ ^ ( I ? Xr,k^ , „. | ) = Z), 

where the average (. . .) is taken over the random choices made only above the (n—k)- 
th stage of the cascade. The factor Pq{k) can be calculated exactly (see below) and 
we propose to use it to compensate for the errors introduced by the reconstruction 
procedure. In particular the reconstructed moments entering (^ will be shifted by : 

^q,reconstructed ^ ^q,reconstTucted ~ log(Pij(^))- (13) 

We will now determine the explicit form of the correction Pq{k). By the definition 
of the a model, the correction Pq{k) can be calculated just by evaluating : 





in the a model modified by taking the starting bin Xq'^ = 1. 

First it is easy to see that for q = 1 there is no correction pi{k) = 1. This is due 
to d^). Also all corrections vanish for k = : 

pM = 1- (15) 

The appearance of a correction for g > 1 comes from the fact that the 'number' of 
particles in this model has a nonzero dispersion. 

Consider first the case of g = 2. We will now split the bins (xj's) appearing in 
(III) into a left half {i < 2''-^) and a right half {i > 2''-^): 



(16) 
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Using the fact that the left and right bins are independent one gets the recurrence 
relation: 

P2ik) = - {paa^ + Pbb^)p2ik - 1) + -. (17) 
z ^ ^ ^ z 

This can be solved together with the initial data (|1^), to yield a closed form solution: 

In general one can obtain the recurrence relation for general q in exactly the 
same way: 

P.(^) = ^ E ( ■ ) d^d,.^ ■ Pri^ - - 1) (19) 



2=0 



where 



5 Discussion 



di = Paa' + Pbb\ (20) 



We have performed numerical simulations of the a-model in order to test the im- 
proved single data analysis in practice. In Fig. la-d the histograms of the corrected 
(with the shift (|19]) taken into account) and standard (without the correction (|19|)) 



values of intermittency exponents ip2, v^a are plotted for 90000 generated cascades 
of 5 and 10 steps. The peaks with the correction included are significantly closer 
to the theoretical value. The dispersion of the distribution estimated directly from 
the observed peak, for the "corrected" histogram is smaller than the dispersion of 
the "standard" one. It decreases with the number of cascade steps. The numerical 
values of "corrected" and "standard" dispersion as a function of the cascade length 
are presented in Tables 1, 2 for 2 different sets of cascade parameters. The corrected 
dispersion is relatively small, and it allows to distinguish between the cascades with 
different parameters (Figs. la-d). 

The influence of the correction (|19D on the value of the intermittency exponents 
obtained from averaging over the ensemble of events ('center of mass' of the his- 
togram) was also investigated. The results are presented in Tables 3,4 for 2 different 
sets of cascade parameters. The estimation of intermittency exponents for the cor- 
rected case is much better than for the standard one. 

In the preceding, the formula for the correction (see e.g. (|18|)) depends on the 
values of the parameters a, b of the a-model. In practice, however, one would like to 
implement some sort of model independent correction. A possible way of doing this 
is to use the fact that the corrections \ogp2{i) and log p3{i) seem to change most 
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dramatically in the first few steps of the reconstruction procedure (near the 'end' of 
the cascade). After that they seem to stabilize at some constant value. This would 
suggest using just the reconstructed moments near the beginning of the cascade in 
the fit (^). In practice, however, this might perhaps suffer from low statistics and 
large fluctuations. 

An alternative procedure would be to first determine the parameters a and h 
using the standard (uncorrected) method, and then substitute those parameters 
into (|T^ and use the improved analysis to obtain a better approximation of the 
exponents. One could repeat this until the result no longer changed. 

6 Conclusions 

Our conclusions can be summarized as follows : 

(a) the value of intermittency exponent estimated from the maximum of "cor- 
rected" histogram moves closer to the theoretical value, 

(b) the dispersion of the distribution estimated directly from the observed peak 
for the "corrected" histogram is smaller than the dispersion of the "standard" one, 

(c) the corrected value of intermittency exponent obtained after averaging over 
the sample of events estimates the theoretical value better than in the standard case, 

(d) a possible procedure of improving the analysis without the knowledge of 
a-model parameters is proposed. 
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Figure captions Histograms of the intermittency exponents 02 (left column) and 
03 (right column) simulated for the set of parameters a = 0.8, b = 1.1 (upper row) 
and a = 0.5, b = 1.5 (lower row) in 90000 events for 5 and 10 cascade steps. The 
wider curves correspond to 5 stages of the cascade. 'Solid' curves represent the 
histograms with the correction (19) taken into account. 
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Table 1. Standard and corrected intermittency exponents (determined from the 
position of the maximum of the histograms) and their dispersions (errors) for 
a — 0.8, h — 1.1 and n — 5, . . . , 10 cascade steps. Theoretical values for inter- 
mittency exponents are (fi2,theor — 2.85 x 10"^ and (ps^theor — 8.13 x lO"^. 



(fi = 10-' X 


5 


6 


7 


8 


9 


10 




1.90 ±0.89 


2.00 ±0.82 


2.26 ±0.75 


2.49 ±0.66 


2.52 ±0.59 


2.46 ±0.52 


^2,corr 


2.66 ±0.75 


2.79 ±0.66 


2.66 ±0.56 


2.85 ± 0.46 


2.85 ± 0.43 


2.85 ±0.36 


^3 


5.66 ±2.46 


5.66 ±2.46 


6.48 ±2.05 


6.64 ± 1.72 


6.81 ± 1.72 


6.72 ± 1.49 


^3,corr 


7.79 ±2.13 


7.63 ± 1.81 


7.62 ± 1.64 


8.00 ± 1.40 


8.12 ± 1.15 


8.12 ±0.98 



Table 2. Intermittency exponents (determined from the position of the maximum of 
the histograms) and their dispersions (errors) for a = 0.5, 6 = 1.5 and n = 5, . . . , 10 
cascade steps. Theoretical values for intermittency exponents are (p2,theor = 3.22 x 
10-1 and ip3,theor = 8.07 x 10-^. 



ipi = 10-1 X 


5 


6 


7 


8 


9 


10 




2.00 ± 1.00 


2.09 ±0.92 


2.43 ±0.74 


2..61±0.70 


2.44 ± 0.74 


2.26 ±0.70 


^2,corr 


3.13 ±0.79 


3.13 ±0.74 


2.96 ±0.61 


2.87 ±0.57 


3.05 ±0.52 


3.00 ±0.48 


^3 


4.52 ±2.23 


5.31 ±2.23 


5.83 ±2.02 


5.90 ± 1.83 


6.03 ± 1.83 


5.70 ± 1.70 


^3,corr 


7.53 ± 1.90 


7.93 ± 1.83 


7.66 ± 1.57 


7.53 ± 1.31 


7.66 ± 1.31 


7.66 ± 1.18 
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Table 3. Standard and corrected intermittency exponents and their dispersions (er- 
rors) for a = 0.8, 6 = 1.1 and n = 5, . . . , 10 cascade steps obtained after averaging 
over the sample of 90000 events.Theoretical values for intermittency exponents are 
(P2,theor = 2.85 X lO"' and (/?3,tw = 8.13 x 10"'. 



(fi = 10-' X 


5 


6 


7 


8 


9 


10 




2.16 ±4.34 


2.45 ± 1.51 


2.57 ±0.74 


2.63 ±0.66 


2.68 ±0.60 


2.72 ±0.54 


^2,corr 


2.90 ±0.70 


2.89 ±0.58 


2.88 ±0.49 


2.88 ±0.42 


2.87 ±0.36 


2.87 ±0.31 


^3 


6.54 ±4.80 


7.03 ±2.62 


7.33 ±2.11 


7.53 ± 1.90 


7.68 ± 1.73 


7.77 ± 1.58 


^3,corr 


8.32 ± 2.00 


8.29 ± 1.67 


8.25 ± 1.40 


8.23 ± 1.18 


8.22 ± 1.02 


8.20 ±0.89 



Table 4. Intermittency exponents and their dispersions (errors) for a = 0.5, b = 1.5 
and n = 5, . . . , 10 cascade steps obtained after averaging over the sample of 90000 
events. Theoretical values for intermittency exponents are (p2,theor = 3.22 x 10^^ 
and (p3,theor = 8.07 x 10"^ 



= 10-1 X 


5 


6 


7 


8 


9 


10 




2.33 ± 1.2 


2.50 ±0.95 


2.62 ±0.82 


2.69 ±0.76 


2.77 ±0.74 


2.81 ±0.72 


V 2,(.'0/7- 


3.20 ± 0.70 


3.17 ± 0.63 


3.15 ± 0.57 


3.15 ± 0.52 


3.11 ±0.17 


3.1 1± 0.13 




5.78 ±2.36 


6.13 ±2.06 


6.38 ± 1.90 


6.55 ± 1.81 


6.71 ± 1.78 


6.81 ± 1.75 


^3,corr 


8.16 ± 1.71 


8.05 ± 1.51 


7.96 ± 1.36 


7.90 ± 1.24 


7.86 ± 1.13 


7.84 ± 1.06 
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Figure 1: a-d 
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